ON THE COMPUTATION OF THE OMEGA INVARIANT OF A NUMERICAL 
SEMIGROUP BY OPTIMIZING OVER AN EFFICIENT INTEGER SET 

VICTOR BLANCO 

Abstract. In this paper we present a mathematical programming fornmlation for the oj invariant of 
a numerical semigroup for each of its minimal generators. The model consists of solving a problem 
!^»J , of optimizing a linear function over the efficient set of a multiobjective linear integer program. We 

fvj ' offer a methodology to solve this problem and we provide some computational experiments to show the 

applicability of the proposed algorithm. 
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1. Introduction 



r ^ , Applying discrete optimization usually means solving real-world problems in industry, transportation, 

^**^ I location, etc, but there are many problems arising in abstract mathematics, specially in commutative 

algebra, that can be formulated as integer programming problems and that have centered the attention 
of many researchers. In particular, problems in classical number theory are specially interesting and only 
a few papers appear in the literature where operational researchers deal with these type of problems. 
In the recent years, there have appeared several methods based on algebraic tools to solve (single 
^SJ , and multi-objective) discrete optimization problems as Grobncr bases [5J|5] or short generating functions 

t^ • [31 El I12[ I13[ 121) . amongst others. In this paper we intend to give an inverse viewpoint, that is, how 

^S| , optimization tools can help to solve algebraic problems and then getting closer both mathematical fields, 

algebra and operational research. 

Here, we study a specific problem in a basic algebraic structure: numerical semigroups. It is a simple 
f~^ ' framework where developing optimization tools to make computations that are usually done by brute 

^— ^ ■ force. 

A numerical semigroup is a set of non-negative integers, closed under addition, containing zero and 
such that its complement in N is finite. Details about the theory of numerical semigroups can be found 
in the recent monograph by Rosales and Garcia- Sanchez [55]. 
Cu I It can be checked that many of the main notions in the theory of semigroups are, in fact, defined or 

characterized as discrete optimization problems. For instance, for a numerical semigroup 5, its multi- 
plicity, m(S'), is the least integer belonging to S, that is 

m(S') = min{s G S}. 
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Note that, if the reader is not famihar with numerical semigroups, may be it may not be possible to see 

why the above problem is a discrete optimization problem, but any numerical semigroup is minimally 

generated by a finite set of integers {rii, . . . , Up} with gcd(ni, . . . , Up) — 1, and then, any element in S 

can be expressed as a linear combination, with non-negative integer coefficients, of this set. Hence, the 

multiplicity of S can be written as 

p 
m(5') = min{no : \ XiUi — nQ,x E 1/^,no E Z+}. 

In this case, it is clear that m(S') = minjrii, . . . , rip}. 

Other well-known index is the Frobenius number. If S* is a numerical semigroup, its Frobenius number, 
F(5), is the largest integer not belonging to S. This index is well defined since by definition NyS' is finite. 
Then, 

F(S') = max{n G Z\5} 

Selmer proved that when a numerical semigroup, S*, is generated by only two elements a and b, then 
F(5) = a6 — (a + &). However, there have not been many successes in finding formulas to compute the 
Frobenius number when the numerical semigroup is generated by more than two integers. Unfortunately, 
this problem cannot be expressed as a linear integer problem. 

The last classical index that we want to mention is the genus of a numerical semigroup. If S is 
a numerical semigroup, the genus of 5, g(S'), is the cardinal of N\5 (its number of holes), that is, 
g(5) = #N\5. 

New arithmetic invariants for commutative semigroups, and then, in particular for numerical semi- 
groups that have recently appeared in the literature are the tame degree, the catenary degree and the cj 
invariant (see [U [SJ [TTJ [T71 [THl HB [HI US] for a complete algebraic description of these indices). These 
arithmetic invariants come from the field of commutative algebra and factorization theory, and their 
definitions are not intuitive from a combinatorial viewpoint. In particular, the uj invariant allows to 
derive crucial finitencss properties in the theory of non- unique factorizations for semigroups (tamcness) . 
However, in [5] the authors give characterizations of some of these concepts that allow the optimization 
approach giving here. 

The goal of this paper is to present a new method for computing the lo invariant of a numerical 
semigroup by optimizing a linear function over the non-dominated solutions of a multiobjective integer 
programming problem. In contrast to usual integer programming problems, in multiobjective problems 
there are several objective functions to be optimized. Usually, it is not possible to optimize all the 
objective functions simultaneously since objective functions induce a partial order over the vectors in the 
feasible region, so a different notion of solution is needed. A feasible vector is said to be a non-dominated 
(or Pareto optimal ) solution if no other feasible vector has component-wise larger objective values. The 
evaluation through the objectives of a non-dominated solution is called efficient solution. 

There are nowadays relatively few exact methods to solve general multiobjective integer and linear 
problems (see |15)). Some of them combine optimality of the returned solutions with adaptability to a 
wide range of problems [31M22| . Other methods, such as Dynamic Programming, are general methods for 
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solving, not very efficiently, general families of optimization problems (see |29]). A different approach, as 
the Two-Phase method (see [lH])j looks for supported solutions (those that can be found as solutions of 
a single-objective problem over the same feasible region but with objective function a linear combination 
of the original objectives) in a first stage and non-supported solutions arc found in a second phase 
using the supported ones. The Two-phase method combines usual single-criteria methods with specific 
multiobjective techniques. 

Optimizing over an efficient set consists then of maximizing (or minimizing) a single objective function 
over the solutions of a multiobjective problem. Since it is not known the structure of the solution set of 
that problem, different techniques that those used to solve single or multiobjective problems are needed 
if we do not want to enumerate all the efficient solutions and then evaluate the function over all of them. 
Although there have appeared some papers analyzing this problem when the problem is continuous (see 
for instance, dj or [30]), only a few study the discrete case (pQ [501 US])- We apply an adaptation of 
the algorithm presented in }20j to compute the a; invariant of a numerical semigroup. We choose that 
methodology since the approaches given in [l] and |23| does not fit the requirement of our problem. 

Although there is no need of computing the w invariant of a numerical semigroup in a short time, 
the hardness of computing this index even for a small number of generators makes it difficult to analyze 
the invariant. Furthermore, algebraic researchers are showing an interest in studying certain combinato- 
rial properties of non- unique factorizations (in particular by the w invariant), but this difficulty makes 
impossible to go further in the study of this index. In [5] the authors give an algorithm to compute 
the uj invariant of a numerical semigroup by using techniques from factorization theory in monoids, and 
they give a short list of small problems that can be solved by using that methodology. Wc provide here 
a new method that avoid the complete enumeration of the non-dominated solutions of the subjacent 
multiobjective problem, and then, we are able to compute more efficiently the lo invariant of a numerical 
semigroup. This allows us to obtain computations that cannot be done up to the moment by using the 
brute force algorithm and then, better analysis of the behavior of this invariant. In the last section of 
the paper wc show that our method, not only improve the CPU computation times, but the sizes of the 
problems that can be solved. However, in this paper we do not underestimate the algorithm in [2]. Our 
algorithm seems to be able to compute the invariant for larger instances, but the approach in [5] allows 
to give interesting conclusions for some special families of numerical semigroups. We center here in the 
computational side of the problem. Other advantage of formulating the omega index as a mathematical 
programming problem is that a better analysis of the constraints and the objective function may improve 
the general running of the method. 

In Section [2] we introduce the preliminary definitions and results for the rest of the paper. There, we 
describe the objects under study, the numerical semigroups, and the index we are interested to compute, 
the w-invariant. An algorithm for optimizing a linear function over an efficient integer set is described 
in Section [3J and it allows to compute efficiently the w-invariant of a numerical semigroup. In Section 2] 
we show some computational tests done to check the efficiency of the algorithm. This paper is almost 
self-contained and written to be accessible for both algebraists and operational researchers. 
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2. Preliminaries 

We begin this section introducing the structure under the study of this paper. A numerical semigroup 

is a subset S' of N (here N denotes the set of non-negative integers) closed under addition, containing zero 

p 

and such that N\S' is finite. We say that {ni, . . . , rip} is a system of generators of S* if 5* = { ^ 

i=i 
Xi e N, i = 1, . . . ,p}. We denote S = (rii, . . . , Up) if {rii, . . . , Up} is a system of generators of S. Note 

that any numerical semigroup has an unique minimal system of generators, in the sense that no proper 

subset of it is a system of generators. Then, we assume through this paper that if 5' = (ni, . . . , Up) then 

{ni, . . . , Up} is the minimal system of generators of S and ni < • • • < Up. The cardinality of the minimal 

system of generators, p, is called the embedding dimension of S. 

Note that, for S ~ {ui, . . . , np) a numerical semigroup we can define a honiomorphism tt : N^ — ;■ 5" as 

p 
7r(a;i,.. .,Xp) = y^^XjUj 

i=l 

usually called the factorization homomorphism. For any element s in 5 we define Z(s) as the set 7r^^(s) = 

p 
{x G N^ : > XiUi = s}. The factorization monoid is the set Z(S') = it^^{S). 

4=1 

For these semigroups we are interested in computing an index that measures how far away generators 
of a numerical semigroup are from primes: the u! invariant (see [THl HZl UHl HI] for further details). 

Definition 2.1. Let S = {ni, . . . ,np) be a numerical semigroup. For s ^ S, let uj{S,s) denote the 
smallest iV G Nq U {oo} with the following property : 

p 
For all n G N and si, . . . , s„ G S", if TJ s^ — s G S, then there exists a subset fl C {1, . . . , n} such 

i=i 
that I f^ I < AT and 

^ Sj - s G S" . 
J en 
Furthermore, we set 

uj{S) = max{a;(5', n^) : i = 1, . . . ,p} G N. 

Note that uj{S) < cx) since any numerical semigroup satisfies the ascending chain condition for w-ideals 
(see [5] and the references therein for further details). 

It is not direct to connect the above definition with optimization theory. However, the following 
result, proved in [5], is crucial to develop the integer programming method proposed in this paper. Here, 
if A C N^, for some p < 1. we denote by Minimals A the set of minimals elements in A with respect to 
the component- wise order in N^. 



Theorem 2.2 ([5]). Let S = (ni, . . . ,rip) be a numerical semigroup. 

p 
1. For every s E S , (jj{S, n) — maxj^^ Xi : x G Minimals Z{n + S)}, 



1=1 



2. uj(S) = max{ > Xi : x E Miniinals(Z(7ii + S)) for some i — 1, . . . ,p\ . 
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P 

Wc denote by |a;| = ^^ Xi the length of x £ N^. From the above result, compute the lj invariant 

consists of maximizing the lengths of the minimal elements in the factorization monoid. 

In the next section we connect the above result with the problem of optimizing a linear function 
over the efficient set of a multiobjective integer linear problem. In what follows, we recall the general 
structure of the optimization problems that are involved in the mathematical programming formulation 
of the omega invariant. This part of the paper can be skipped if the reader is familiar with this theory. 

An integer linear programming problem consists of finding the largest (or smallest) value of a linear 
operator (objective function) c : Z" — )■ Z when c only takes values over a feasible region defined by a 
system of linear equations: 



aiixi + ai2X2+ 



+ ainXn = bi 



amlXl +ajn2X2 + 



where a.^ ,bj E Q for all i G { 1 , 



I ^mnXj] 



i} and j G {!,..., m}. If we denote hy A = {aij) G 



b ~ {bj) G Q"\ this problem can be written as 



and by 



(IP) 



max c(x) 
s.t. 



Ax =6 

X G Z" 



It is well-known |22] that the complexity of solving the above problem is NP-complete, and difficult 
to solve even for a small number of variables (see [S]). However, when the constraints {Ax ~ b) have 
a concrete structure, or the formulation comes from a problem where more information is known, it is 
possible to develop exact (or heuristic with high precision) algorithms that allow to solve efficiently these 
problems (see |24j ) . Recently, some papers have appeared where algebraic tools are used to solve problem 
(Up]) (see lilllllT]). 

If we consider in (jIP[) a linear operator C : Z" — > Z'^ (with k > 2) we get a integer linear multiobjective 
optimization problem: 



(MIP) 



max( or min) C{x) 
s.t. 



Ax 

X eZ'i 



Since, in this case, objective function defines a non-total order over the set of solutions of Ax = b with 
X G Z", it is needed to define what we understand about solving (jMIP|) . We say that a vector x G Z" 
with Ax = 5 is a non- dominated (or Pareto optimal) solution if there is no other vector y G Z" with 
Ay = b and such that C{y) < C{x) and C(a;) / C'(y) where < denotes the component-wise order in 
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ll^ . Then, the solution of (|MIPp is a set of feasible solutions that are not comparable and such that we 
cannot find better solutions with respect to the partial ordering induced by C. In the case we want to 
minimize instead of maximize, the ordering is adapted conveniently. 

In general, the complexity of solving multiobjectivc integer problems is (tP-hard and just a few algo- 
rithms appear in the literature to solve exactly this kind of problems. The recent works by Blanco and 
Puerto [B], De Loera et. al [13] and Blanco and Puerto [7] present algebraic approaches for solving these 
problems. An interesting monograph for multiobjective optimization is |15| . 

In the next and last step we recall the goal of optimizing over an efficient integer region. Note that 
given the set of non-dominated solutions of (jMIP[) we may be interested in finding those non-dominated 
solutions that maximize a linear operator. Then, we arc interested in solving the following optimization 
problem: 

(OES) "^"" '^'^ ^ , 

s.t. a; is a non-dominated solution of (|MIPp 

It is not possible, in general, to express the constraint of the above problem as a linear system of equations. 
Actually, as far as we know, there are no descriptions for the set of solutions of a multiobjective problem, 
and then, to solve (|OES|) one may think of solving first the multiobjective problem, and then, selecting 
those solutions with largest values of c. This process is not effective and it makes it difficult when the 
efficient set has a large size. 

However, there exists a few algorithms to solve (jOESp without enumerating the complete set of non- 
dominated solutions of (jMIPp . In particular, we will adapt the algorithm presented in [20] to solve the 
problem of optimizing over the efficient integer set related to the w-invariant. The following result states 
the structure of the optimization problem that solves the w-invariant. 

Theorem 2.3. Let S = {ni, . . . , Up) be a numerical semigroup. Then, for each j G {1, . . . ,p}, Uj(S, Uj) 
is the solution of the following OES problem: 



(OES°) 



X € min C{x) = {xi, . . . ,Xp) 



max c[x} = 2^. X 

i=l 

s.t. 



" s.t. 



p p 

i=l 1=1 

x,y ez\ 



Proof. Observe that the elements in Minimals(Z(r7,j +5)) arc those minimal elements with respect to the 

p 
component wise ordering in Z(r7,j + S*), and Z(r7,j + 5) = {{xi, . . . , Xp) :y UiXi £ Uj + s, for some s G 



i=l 



S} = {(xi, . . . , Xp) : 2, "-i Xi — Uj + N^ Ui Hi for some yi G N}. D 
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Problem ( OES" I consists of maximizing the lengths of the non-dominated solutions of a multiobjective 



problem. We denote by (SMIPj) that subjacent multiobjective problem 

From the formulation of the w-invariant given in the above result, it is clear that if we denote by 
Cj; the p-tuple having a 1 as its jth entry and zeros otherwise, for each j S {1, . . . ,p}, Cj G Z(?ij + S) 
since we can always write rij ~ y ■ rii + 1 ■ rij and (0, . . . , 0) is not a feasible solution. Consequently, 

if the multiobjective subproblem has more non-dominated solutions they must have Xj = to be not 
comparable with Cj . The next lemma states that these unit vector cannot be solution of the multiobjective 
subjacent problem. 

Lemma 2.4. Let S — (ni, . . . , rip) be a numerical semigroup. Then, uj{S, n.j) > 1 for all j G {1, . . . ,p\. 

Proof. From Definition 12.11 ijj{S,nj) > 1. Let us see that uj{S,nj) ^ 1. Assume that uj{S,nj) ~ 1, then, 

71 n 

if 2^, Si — Uj G S, there exists some Sk such that Sk — rij G S. That is, if Uj ^ N, ^i + s for some s G 5, 

then Tij = Sk + s' for some s' G S, for any si, . . . , s„ and any n G N. It is clear that it is not possible 
since that assumption is the same as stating that the generators of S are prime, but there are no primes 
in a numerical semigroup (see [H]). D 



Furthermore, note that the feasible region in (OES"l is not bounded. However, by solving linear 



integer problems we can obtain upper bounds for each of the variables Xi involved in the problem, 
making bounded the feasible region. For each j G {1, . . . ,p} and k ^ j, let ubjk the optimal value of the 
following integer program: 



s.t. 

p 
(UBjfc) XjHj ~ ^ jjin^ = Uk 

1=1 
2/fc=0 

Xj G Z_|_ , y G Z^ 



The solutions of ( UBjfc I are those feasible solutions of (|OES|) in the form Xj Cj . An upper bound for each 
variable in (OES. ) is then ubj = ma.Xk{ubjk}. 

Note that if a; is a non-dominated solution of the multiobjective subproblem in (jOESp . and Xj ^ 0, 
then Uj = 0, and also, if yj ^ 0, then Xj = 0. It is clear since if Xj and yj have nonzero values, then, 
we can always find a feasible solution [x' ,y') where cc'- = max{0,a;j — j/j}, y' = max{0,yj — .tj}, and 
x[ = Xi and y[ = yi for i ^ j. This result is also proved in [S] from an algebraic viewpoint, in terms of 
the supports of x and y. 

As a direct consequence of Lemma 12.41 Theorem 12.31 and the comment above we get the following 
result: 
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Corollary 2.5. Let S = (711, . . . , 7ip) be a numerical semigroup. Then, for each j G {1, . . . ,p}, Uj{S, rij) 
is the solution of the following OES problem: 

n 

max y Xi 

S.t. 



(OES,) 



X G min(xi, . . . ,Xp) 

s.t. 

p p 

y^ XiUi - ^ yiUi = n-j 

i=l i=l 

Xi < ubi 

x,^0 

x,y £Z^ 

Now, computing the omega invariant consists of solving a problem of optimizing a linear function over 
the efficient region of a bounded integer linear multiobjcctive problem. In the next section we describe a 



methodology for solving ( OESj ) 



3. An algorithm for computing uj[S,nj 



In this section we adapt the procedure given in |20j for solving (OESj), and consequently compute 



the omega invariant of a numerical semigroup. For the sake of that, we first give some previous results 
related to each of the steps of the algorithm. 

Let j G {1, . . . ,p}. The first step of the proposed procedure consists of stating that we effectively need 



to design an algorithm to compute uj{S,nj). With this aim, we consider a relaxed problem of (OESj) 



consisting of obviating that the feasible solutions arc just the minimal non-dominated solutions of the 
multiobjcctive problem, but the feasible solutions of (SMIPj). This problem is the following 



(R.) 



max N Xi 




i=i 




s.t. 






p 




y ] xiHi 




i=l 




Xi < ub 




Xj = 



p 

i=l 



,yeZP. 



If we solve (R, ) and the optimal solution is a non-dominated solution of (SMIP,), then, that solution is 



the solution of (OESj ). Furthermore, if (R, ) is infeasible, (OESj ) has no solution. The next result gives 
us some information about the above problem. 



Lemma 3.1. Problem (Rj ) is feasible. Furthermore, the optimal solutions of (R, ) are not non- dominated 
solution o/(SMIPj). 
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Proof. By the definition of tlie a; invariant and tlie cliaracterization given in Theorem 12. 3[ it is clear that 



(Rj I cannot be infeasible. Moreover, the problem is bounded due to the finiteness of a; (5*) (see comment 



after Definition 12.1 



Let (x* , y*) an optimal solution of (Rj ). Assume first that x* = ubi > for all i = 1, . . . ,p with i j^ j 
{x* = 0). Then, by the comment done before Corollary 12.51 y* — for all i ^ j, and then, by feasibility 
Si=^7 -^i ^i ~ ^j' that is not possible since we are supposing that {ni, . . . , Up} is a minimal system of 
generators of S. Then, there exists at least one k ^ j such that x^ < ubk- 

Assume now that (a;*, y*) is a non-dominated solution of (SMIPj), and let x = x* + {ubk — x^.) e^ and 
y = y* + {ubk — 2:^)efc. By definition {x,y) is a feasible solution of (SMIPj) and x > x* , and consequently 
J2^=i ^ > 12^=1 ^* contradicting the optimality oi x* . D 



The above result justifies the development of an algorithm that allows to improve the optimal solution 
of (Rj ) to obtain the optimal value of (UESJ). In the next step we search for feasible solutions of (SMIPj) 



that are non-dominated and that dominate the optimal solution of (Rj I. For that, we use the following 
adaptation of the result by Ecker and Kouada [ll] . 



Lemma 3.2. Let x* be an optimal solution of (Rj) and {s^x,y) an optimal solution of the following 
problem 



n 

max 2, Xi 

i=l 

s.t. 



(EK,(a:*)) 



Xi + Si = x* i = 1, . . . ,p 

p p 

y^ XiUi - ^ yiUi = Uj 

i=l i=l 

Xi < ubi i = 1, . . . ,|) 

Xj ==0 

x,y e TL\. 

Then, x is a non- dominated solution o/(SMIPj) that dominates x* . 

With this lemma, we are able to compute non-dominated solutions of (SMIPj) from the solutions of 
(Kj\. However, these solutions may not be the ones that maximize the lengths in (jOESj ). 



To move through feasible solutions of our problem to obtain the wished solutions, we use the idea 
presented in |20) . where the author gives a procedure to generate new non-dominated solutions. This 
procedure uses a result given in the book by Nemhauser and Wolsey [23] . 
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Lemma 3.3. Let x be a non- dominated solution (SMIPj) and {x,y) an optimal solution of the following 
problem 

n 

max \^ Xi 

i=l 
S.t. 



(NW,(x)) 



Xi < Zi{xi - 1) - Mi{zi - 1) i = 1, 



i=i 

Xi < ubi 


j=i 


rij 


Xj =0 
p 






X, 2/ G Z^ 


ze{0,i}p 





i= l,...,_p 



id where Mj = nia.x{xj : UjXj = > n.^ yi, Xi < ubi, * = Ij • • • iP? 2;, y G Z+}- 
Then, x is a feasible solution of (SMIPj) that is not dominated by x. 



By using the above lemma, in case the solution x is not a solution of ( OESj ) , we can generate other 



feasible solution, non comparable with it and such that we can check if it is the one that we are looking 
for or, if it is not, continue iterating to reach it. 



Therefore, if the solution x obtained is not optimal for (OES, ) and we find x by using Lennrra 13.31 x 



can be replaced by x in the procedure of Lemma l3.2l to obtain a non-dominated solution that dominates 



X. If this new non-dominated solution is not the optimal solution of ( OESj ) we use again Lemma 



to 



obtain a non comparable solution with this. But it may occur that this process cycle obtaining x again. 
However, this can be avoided by applying the following extension of Lemma 13.31 and also the reader can 
find it in [53], that allows to obtain a feasible solution of (SMIPj) that is not comparable with a given 
list of non-dominated solutions x^, . . . ,x'^ . 
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Lemma 3.4. Letx^, 
following problem 



, , x^ be non-dominated solutions o/(SMIPj) and {x,y) an optimal solution of the 



max y Xi 



s.t. 



x,<z^ixf^l)^M,{z^-l) 



p 



(NW,(si,...,x^)) 



i=l 



XiTll 



P 



Xi < ubi 

Xj =0 

p 

i=l 



i,...,p,fc = 1, 
i,...,p 

: 1,...,S 
: 1 S 



id where Mj = niax{xj 



ij^j 



p 

E 



niyi.Xi < uhi.i = 1, . . . ,p,x,y e Z^}. 



Then, x is a feasible solution o/(SMIPj) that is not dominated by x^, 



,x . 



With this method we are able to find a feasible solution that is not comparable with any of the solutions 
found in the previous iterations (those that we know that are not optimal for (jOESp ). We can repeat 
the process until we get the desired solution, but when do we know that wc have got the solution of our 
problem? To ask this question we analyze the last step of the described procedure. 



Assume we have obtained a non-dominated solution of SMIPj, x. If we solve (NWj(a; , . . . , a;^) ) and 



we get X as optimal solution with X^iLi ^i *^ X^iLi ^i then, x is the solution to our problem since among 
all the possible non-dominated solutions, the one that has maximum length has smaller length than x. 



Moreover, if (NWj(a:;^, . . . , x") ) is infeasible, clearly x is our solution. 



The general procedure is describe in the pseudo-code of Algorithm [T] 
Algorithm 1: Computing the omega invariant by optimizing over an efficient set. 



Input 



A minimal generating system of a numerical semigroup S: {ni, 
je {l,...,p}. 



. . , lip 



} and 



Solve ^■. x^. 

while Lu{S,nj) is not obtained do 



• Set X* = x'- and solve (EKj(x*)): x' 

• Solve (|NWj(r,...,x'*)|: x'. 



if (NWj(x^, . . . , x^) I is infeasible or Yl^i=i ^i < X]f=i ^\ then 

I x' = uj{S,nj) 
else 

L Set x'+^ = X and I -.^ I + 1 



Output 



uj{S,nj 



A direct consequence of all the above comments about the proposed procedure is the following. 
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Theorem 3.5. Let S ~ {ni, . . . ,np{ be a numerical semigroup. Algorithm{l\computes uj{S,nj) in a finite 
number of steps. As uj{S) = max{a;(5', n^) : j g l,...,p}, Algorithm[l\ also allows the computation of 

The following example illustrates the usage of Algorithm [1] 

Example 3.6. Let S = (6, 13, 14). The corresponding bounds are ubi = 9, w&2 = 2 and ub^ = 4 and 
the constants for (NWj(x^, . . . , a;*) ) are Mi = 9, M2 = 8 and M3 = 9. The optimal solution for Ri is 
x^ = (0, 9, 9). Solving now (EK(a;-'^)) we obtain x^ = (0, 2, 0), and finally, the optimal solution of NW(2;^) 
is X = (0, 1, 9). Since X]i=i a^i = 10 > J2i=i ^i ~ 2, we continue with the algorithm with x^ = (0, 1,9). 

Computing now the solution of (EK(a;^)) wc get x'^ — (0, 0, 3) and 'NW{x^,x'^) is infeasible, so uj{S, 6) = 
3. 

Analogously we obtain uj{S, 13) = 9 and ^(5", 14) = 7, being uj{S) = max{3, 7, 9} = 9. 



Note that further improvements can be done in the above algorithm by incorporating bounds for the 

p 

w invariant at each step of it. For instance, we can initially set uq ~ Q and vq ~ y x^, the lower and 

upper bounds for uj{S,nj), respectively. Then, update uq with the optimal value of ( |EKj(a:*) ) and vq 
with the optimal value of (NWj(a;^, . . . ,x*) I when they are improved. With this, we have that either 
when uo = uo or the optimal value of (NWj(x^, . . . ,5;^) I is smaller than uq, the solution is the last 
solution obtained by solving (EKj(a;*) ). We have included these bounds to our implementation for the 
computational results, although some others could be added by incorporating these bounds conveniently 
to the constraints of the problems. 

Moreover, if the Apery sets of the numerical semigroup are known, the following lemma, that appears 
in [5], helps us to simplify the problem of searching minimal solutions. 

Lemma 3.7 ([5]). Let S ~ {ni, . . . , rip) be a numerical semigroup, j G {1, . . . ,p}, and (x*, y*) G (SMIPj). 
p 

Then, if x* ^ 0, ^ n^ y^ G Ap{S, ni). 

fc=i 

p p 



Proof. Choose i such that x* ^ 0. Since x* is feasible for (SMIPj), \^ n^- a;^ ~ Uj + y^ n^yl,. From 

fc=i fc=i 

p 

the minimality of (.t*, y*), x* — Bj cannot be minimal, that is, > n^xl. — n,; cannot be expressed as a 

fc=i 
p p 

integer linear combination of rii, . . . , Up, or equivalently, 2, '^fc ^t ~ "-j ~ "-i = ^j + / "-fc Vk ~ ^j ^ i^- 

fc=l k=l 

p 
This implies that V^ Uk yl G Ap(5, nt). D 



fc=i 



Then , with the above result, the equation defining the feasible region of (SMIPj) can be written as 

p 
(3.1) ^ Uk Xk = Uj +w w e[J Ap{S, Hi) 

fc=i i^j 
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Since the above constraint is not linear, we can add part of the information contained on it by incor- 
porating inequalities in the form J2'i=i ^i Vi — '^^^Ur=i Ap('5'i ^i) or X^f^i "-* Vi — ™i^Ur=i ^P('S'j «i) = 
niin{ni, . . . , Up} 

4. Computational experiments 

In this section we present the results of some computational experiments done to analyze the applica- 
bility of the algorithm proposed in Section |31 For the sake of that we have generated a set of instances of 
numerical semigroups with different numbers of minimal generators. From the computational viewpoint, 
up to the moment, we can only compare our algorithm with a recent implementation done in GAP |10j 
for [5] (available upon request), and that consists of enumerating all the feasible solutions of (SMIPj), 
comparing them to determine which are the minimal solutions, and finally selecting those with largest 
length. 

Our algorithm has been implemented in XPRESS-Mosel 7.0 that allows to solve the single-objective 
integer problems involved in the resolution of the lo invariant, by using a branch-and-bound method and 
nesting models by calling the library mmjobs. The algorithms have been executed on a PC with an Intel 
Core 2 Quad processor at 2x 2.50 Ghz and 4 GB of RAM. 

It is clear that the complexity of the algorithm depends of the dimension of the space (embedding 
dimension of the numerical semigroups) and the size of the coefficient of the constraints. Then, we have 
generated for each p G {5, 10, 15, 20}, five instances of p non-negative integers ranging in [2, 1000]. Note 
that larger instances are difficult to find since even for p = 20, selecting 20 no n- negative integers with 
great common divisor equal to one and such that they are a minimal system of generators of a numerical 
semigroup is not an easy task. We have used recursively the function RandomListForNS of GAP until 
we found the list of integer with the above requirements. In Table [T] we show the battery of numerical 
semigroup for which we have run the algorithms. In the first column we write the name we give to the 
semigroup, in the second column the set of minimal generating system of the numerical semigroup, and 
in the last column the embedding dimension of the semigroup. 



name 


minimal generating system 


embdim 


S5(l) 
S5(2) 
S5(3) 
S5(4) 
S5(5) 


{20,354,402,417,429} 
{7,292,359,645,755} 

{5,86,99,148, 152} 
{41,65,155,317,377} 
{28,55,125,233,590} 


5 


SlO(l) 
S10(2) 
S10(3) 
S10(4) 
S10(5) 


{43, 63, 68, 108, 120, 135, 142, 150, 177, 224} 

{15, 46, 58, 89, 108, 114, 117, 126, 130, 173} 

{20, 22, 24, 26, 54, 77, 83, 89, 93, 95} 

{96, 131, 136, 171, 173, 239, 278, 287, 364, 483} 

{146, 173, 207, 359, 426, 548, 604, 606, 657, 702} 


10 
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S15(l) 
S15(2) 
S15(3) 
S15(4) 
S15(5) 


{36, 47, 65, 79, 82, 84, 91, 96, 100, 109, 121, 124, 134, 139, 169} 

{46, 115, 155, 286, 289, 341, 342, 348, 393, 436, 445, 449, 504, 527, 584} 

(40, 84, 126, 130, 132, 135, 142, 152, 165, 183, 217, 221, 229, 273, 323} 

{75, 104, 114, 128, 216, 219, 241, 271, 309, 310, 321, 327, 340, 352, 371} 

{29, 50, 95, 96, 99, 109, 110, 119, 131, 134, 135, 152, 162, 180, 201} 


15 


S20(l) 
S20(2) 
S20(3) 
S20(4) 
S20(5) 


{131, 145, 249, 257, 260, 319, 354, 459, 465, 469, 487, 572, 575, 587, 606, 607, 652, 674, 694, 762} 
{57, 105, 182, 186, 201, 204, 254, 259, 263, 274, 275, 294, 295, 298, 307, 338, 367, 393, 417, 431} 
{85, 298, 333, 342, 349, 358, 401, 415, 462, 480, 556, 569, 583, 609, 619, 708, 710, 752, 821, 853} 
{81, 107, 168, 194, 230, 236, 274, 277, 286, 290, 305, 310, 348, 351, 366, 379, 396, 416, 521, 583} 

{101, 141, 279, 314, 329, 369, 399, 425, 438, 447, 477, 501, 534, 536, 555, 574, 620, 727, 786, 871} 


20 



Table 1: Battery of numerical semigroups for the computational 
experiments. 



In tables [5] and |3] we summarize the obtained results. In the first column we write the numerical 
semigroup (S). The second column is the generator with respect to the w invariant is computed (n^). 
The value of Lo{S,nj) is shown in the third column (in this same column we mark in bold face, for 
each numerical semigroup, the maximum value that coincides with a;(S)). One of the minimal vectors 
that reach the maximum length is shown in the fourth column. The next column, the fifth, shows the 
number of iterations of our algorithm (it). The CPU time (in seconds) needed to compute uj{S,nj) 
with our algorithm (timej) and with GAP (GAPtime - * indicates that GAP was not able to solve the 
problem in less than 2 hours) is shown in columns sixth and seventh, respectively. The eigth shows the 
total CPU time needed to compute ljj{S) (totaltime). The average CPU time to compute each of the 
uj(S, Hj) (avtime) is written in column nineth. The last column is the number of minimal (nondominated) 
elements of the corresponding problem that was computed with GAP (* indicates that GAP was not able 
to compute this number in less than two hours). 

For the computations of the experiments for p = 15, 20 we add to the constraints of the 
multiobjective problem the inequalities J2^=i "■« Vi — ^^^^Uf^i API'S"; "•«) or 
Z^Li "^i Vi ^ minULi Ap(S', Ui) = min{7ii, . . . , rip}. 



s 


Hj 


a;(S,n,) 


min 


it 


timej 


GAPtime 


tottime 


avtime 


#miii 




20 


4 


[0,0,0,0,4] 


9 


0.54 


6.03 






12 




354 


60 


[60,0,0,0,0] 


12 


0.95 


11.35 






14 


S5(l) 


402 


63 


[63,0,0,0,0] 


16 


1.439 


12.54 


5.921 


1.184 


17 




417 


60 


[60,0,0,0,0] 


15 


1.43 


12.68 






16 




429 


60 


[60,0,0,0,0] 


17 


1.55 


12.43 






20 




7 


3 


[0,3,0,0,0] 


10 


0.55 


12.48 






11 




292 


93 


[93,0,0,0,0] 


9 


0.37 


23.72 






11 


S5(2) 


359 


93 


[93,0,0,0,0] 


11 


0.43 


27.33 


2.84 


0.56 


13 




645 


200 


[200,0,0,0,0] 


14 


0.67 


45.92 






15 




755 


200 


[200,0,0,0,0] 


18 


0.81 


75.59 






19 
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5 


2 


[0,0,0,2,0] 


8 


0.285 


1.201 






11 




86 


37 


[37,0,0,0,0] 


11 


0.294 


2.527 






12 


S5(3) 


99 


37 


[37,0,0,0,0] 


11 


0.34 


2.82 


1.69 


0.33 


12 




148 


60 


[60,0,0,0,0] 


12 


0.37 


4.1 






13 




152 


60 


[60,0,0,0,0] 


12 


0.39 


2.29 






13 




41 


14 


[0,14,0,0,0] 


12 


0.893 


5.64 






14 




65 


22 


[22,0,0,0,0] 


13 


0.988 


6.02 






14 


S5(4) 


155 


24 


[24,0,0,0,0] 


16 


1.1 


8.22 


7.39 


1.47 


18 




317 


22 


[21,0,1,0,0] 


22 


2.916 


13.96 






28 




377 


31 


[31,0,0,0,0] 


18 


1.49 


18.7 






35 




28 


10 


[0,10,0,0,0] 


11 


0.5 


10.71 






12 




55 


25 


[25,0,0,0,0] 


8 


0.381 


11.45 






12 


S5(5) 


125 


27 


[27,0,0,0,0] 


13 


0.71 


20.18 


4.719 


0.94 


15 




233 


26 


[26,0,0,0,0] 


13 


0.732 


42.37 






17 




590 


30 


[24,5,0,1,0] 


23 


2.38 


109.38 






48 




43 


5 


[0,0,5,0,0,0,0,0,0,0] 


49 


3.36 


5.41 






58 




63 


8 


[8,0,0,0,0,0,0,0,0,0] 


48 


2.7 


8.61 






65 




68 


8 


[8,0,0,0,0,0,0,0,0,0] 


49 


3.16 


13.18 






69 




108 


7 


[5,0,2,0,0,0,0,0,0,0] 


52 


4.15 


18.26 






81 




120 


8 


[8,0,0,0,0,0,0,0,0,0] 


57 


4.24 


12.65 






94 


SlO(l) 


135 


9 


[9,0,0,0,0,0,0,0,0,0] 


68 


5.95 


15.5 


67.89 


6.78 


108 




142 


9 


[9,0,0,0,0,0,0,0,0,0] 


88 


9.24 


16.75 






125 




150 


7 


[4,2,1,0,0,0,0,0,0,0] 


66 


8.07 


19.85 






116 




177 


9 


[7,2,0,0,0,0,0,0,0,0] 


70 


6.88 


49.26 






149 




224 


9 


[7,0,2,0,0,0,0,0,0,0] 


113 


20.1 


65.16 






246 




15 


3 


[0,0,3,0,0,0,0,0,0,0] 


36 


1.801 


6.64 






45 




46 


9 


[9,0,0,0,0,0,0,0,0,0] 


39 


1.66 


10.32 






48 




58 


10 


[10,0,0,0,0,0,0,0,0,0] 


38 


1.69 


12.23 






50 




89 


9 


[7,0,1,0,0,0,0,1,0,0] 


47 


2.681 


17.33 






68 




108 


15 


[15,0,0,0,0,0,0,0,0,0] 


63 


3.278 


24 






83 


S10(2) 


114 


16 


[16,0,0,0,0,0,0,0,0,0] 


57 


3.07 


28.81 


35.87 


3.58 


78 




117 


15 


[15,0,0,0,0,0,0,0,0,0] 


63 


4.316 


21.65 






88 




126 


16 


[16,0,0,0,0,0,0,0,0,0] 


73 


4.243 


22.48 






99 




130 


22 


[22,0,0,0,0,0,0,0,0,0] 


64 


3.399 


38.59 






98 




173 


23 


[23,0,0,0,0,0,0,0,0,0] 


107 


9.73 


80.49 






161 




20 


4 


[0,0,0,4,0,0,0,0,0,0] 


39 


1.48 


5.1 






43 




22 


5 


[5,0,0,0,0,0,0,0,0,0] 


43 


1.59 


5.41 






45 




24 


5 


[5,0,0,0,0,0,0,0,0,0] 


36 


1.3 


5.41 






45 




26 


5 


[3,2,0,0,0,0,0,0,0,0] 


33 


1.64 


3.49 






44 


S10(3) 


54 


6 


[0,6,0,0,0,0,0,0,0,0] 


52 


2.77 


14.05 


99.49 


9.94 


88 




77 


9 


[7,0,2,0,0,0,0,0,0,0] 


93 


13.27 


26.72 






176 




83 


9 


[6,0,2,1,0,0,0,0,0,0] 


109 


19.41 


33.83 






198 
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89 


10 


[10,0,0,0,0,0,0,0,0,0] 


100 


13.85 


41.18 






219 




93 


10 


[10,0,0,0,0,0,0,0,0,0] 


109 


21.75 


46.17 






254 




95 


10 


[10,0,0,0,0,0,0,0,0,0] 


114 


22.4 


52.46 






251 




131 


7 


[5,0,0,0,2,0,0,0,0,0] 


63 


8.34 


61.44 






102 




136 


6 


[3,1,0,0,2,0,0,0,0,0] 


47 


7.23 


54.38 






88 




171 


6 


[2,2,0,0,2,0,0,0,0,0] 


65 


11.18 


56.92 






102 




173 


7 


[3,1,3,0,0,0,0,0,0,0] 


60 


9.87 


116.22 






118 


S10(4) 


239 


8 


[5,2,0,0,0,0,0,1,0,0] 


83 


16.81 


104.66 


225.55 


22.55 


155 




278 


10 


[10,0,0,0,0,0,0,0,0,0] 


80 


14.93 


129.1 






208 




287 


10 


[10,0,0,0,0,0,0,0,0,0] 


62 


11.628 


128.1 






178 




364 


10 


[7,3,0,0,0,0,0,0,0,0] 


128 


34.053 


227.12 






260 




483 


11 


[9,1,0,0,0,0,1,0,0,0] 


204 


105.146 


497 






427 




146 


8 


[0,6,2,0,0,0,0,0,0,0] 


42 


8.048 


100.82 






70 




173 


10 


[10,0,0,0,0,0,0,0,0,0] 


71 


15.43 


115.39 






99 




207 


10 


[10,0,0,0,0,0,0,0,0,0] 


60 


11.77 


138.87 






82 




359 


12 


[7,5,0,0,0,0,0,0,0,0] 


60 


14.69 


198.246 






152 


S10(5) 


426 


12 


[12,0,0,0,0,0,0,0,0,0] 


77 


16.23 


290.08 


315.14 


31.51 


130 




548 


12 


[0,12,0,0,0,0,0,0,0,0] 


105 


38.525 


470.76 






209 




604 


15 


[15,0,0,0,0,0,0,0,0,0] 


124 


43.81 


499.9 






244 




606 


13 


[13,0,0,0,0,0,0,0,0,0] 


98 


28.4 


422.96 






243 




657 


12 


[0,8,4,0,0,0,0,0,0,0] 


105 


65.01 


558.71 






244 




702 


14 


[14,0,0,0,0,0,0,0,0,0] 


159 


73.19 


718.58 






362 



Table 2: Results of computational experiments for p = 5, 10. 



s 


Hj 


a;(S,nj) 


min 


it 


timej 


GAPtime 


tottime 


avtime 


#min 




47 


6 


[6,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


114 


8.386 


13.78 






129 




65 


5 


[4,0,0,1,0,0,0,0,0,0,0,0,0,0,0] 


112 


10.358 


35.64 






159 




79 


5 


[3,1,0,0,1,0,0,0,0,0,0,0,0,0,0] 


105 


9.392 


7.21 






165 




82 


6 


[0,6,0,0,0,0,0,0,0,0,0,0,0,0,0] 


141 


13.664 


17.93 






184 




84 


6 


[4,1,0,0,1,0,0,0,0,0,0,0,0,0,0] 


112 


11.156 


28.24 






192 




91 


7 


[7,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


101 


6.863 


9.34 






173 


S15(l) 


96 


7 


[4,3,0,0,0,0,0,0,0,0,0,0,0,0,0] 


104 


11.98 


52.225 


612.099 


40.8066 


250 




100 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


187 


26.425 


29.725 






251 




109 


7 


[7,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


129 


12.659 


35.725 






245 




121 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


154 


19.271 


48.225 






307 




124 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


214 


37.796 


203.8 






364 




134 


7 


[5,1,1,0,0,0,0,0,0,0,0,0,0,0,0] 


168 


29.652 


241.9 






383 




139 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


183 


30.122 


199.05 






394 




169 


8 


[5,2,0,0,0,1,0,0,0,0,0,0,0,0,0] 


285 


38.405 


164.625 






680 




46 


5 


[0,1,3,0,1,0,0,0,0,0,0,0,0,0,0] 


83 


8.383 


79.4 






98 




115 


6 


[2,0,3,0,1,0,0,0,0,0,0,0,0,0,0] 


94 


10.454 


112.65 






109 




155 


17 


[17,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


123 


14.728 


139.425 






151 



THE OMEGA INVARIANT BY OPTIMIZATION OVER AN EFF. SET 



17 





286 


15 


[12,0,3,0,0,0,0,0,0,0,0,0,0,0,0] 


137 


20.015 


291.65 






206 




289 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


109 


14.545 


293.575 






190 




341 


17 


[17,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


174 


65.975 


401 






252 




342 


15 


[14,0,0,0,1,0,0,0,0,0,0,0,0,0,0] 


192 


32.986 


406.775 






265 


S15(2) 


348 


15 


[13,0,2,0,0,0,0,0,0,0,0,0,0,0,0] 


193 


113.383 


427.3 


1683.63 


112.242 


291 




393 


20 


[20,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


228 


135.869 


550.35 






320 




436 


25 


[25,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


273 


74.036 


736.575 






413 




445 


24 


[24,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


311 


96.198 


784.625 






434 




449 


19 


[19,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


294 


82.945 


795.45 






425 




504 


22 


[22,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


354 


177.154 


1161.45 






594 




527 


20 


[20,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


367 


166.69 


1345.275 






610 




584 


26 


[26,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


438 


670.269 


1988.875 






737 




40 


3 


[0,1,0,0,1,0,0,0,0,1,0,0,0,0,0] 


89 


7.53 


24.85 






113 




84 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


119 


11.648 


39.95 






139 




126 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


114 


13.082 


67.375 






195 




130 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


117 


12.165 


68.425 






188 




132 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


126 


12.849 


69.1 






186 




135 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


149 


23.85 


72.4 






212 




142 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


156 


20.35 


78.35 






216 


S15(3) 


152 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


176 


22.83 


90.65 


1407.684 


93.8456 


267 




165 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


183 


31.199 


109.15 






305 




183 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


189 


38.537 


134.025 






324 




217 


11 


[9,1,0,0,1,0,0,0,0,0,0,0,0,0,0] 


244 


77.582 


219.075 






446 




221 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


250 


70.445 


229.625 






454 




229 


13 


[13,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


274 


89.679 


253.775 






478 




273 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


353 


133.473 


475.725 






701 




323 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


578 


842.465 


3487.55 






1076 




75 


5 


[0,4,0,1,0,0,0,0,0,0,0,0,0,0,0] 


92 


8.229 


42.65 






123 




104 


6 


[4,0,0,1,1,0,0,0,0,0,0,0,0,0,0] 


96 


10.708 


52.2 






140 




114 


7 


[5,0,0,2,0,0,0,0,0,0,0,0,0,0,0] 


86 


8.044 


54.875 






144 




128 


7 


[5,0,2,0,0,0,0,0,0,0,0,0,0,0,0] 


104 


10.528 


61.075 






147 




216 


8 


[3,0,5,0,0,0,0,0,0,0,0,0,0,0,0] 


151 


23.513 


126.85 






254 




219 


9 


[8,1,0,0,0,0,0,0,0,0,0,0,0,0,0] 


145 


19.017 


126.95 






251 




241 


9 


[9,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


139 


22.934 


156.65 






302 


S15(4) 


271 


8 


[4,3,1,0,0,0,0,0,0,0,0,0,0,0,0] 


188 


52.956 


202.2 


892.038 


59.4692 


357 




309 


9 


[6,1,1,1,0,0,0,0,0,0,0,0,0,0,0] 


211 


62.039 


284.25 






439 




310 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


188 


34.511 


281.525 






429 




321 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


269 


67.673 


320.65 






523 




327 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


185 


43.396 


341.225 






502 




340 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


212 


42.852 


374.325 






509 




352 


10 


[7,3,0,0,0,0,0,0,0,0,0,0,0,0,0] 


262 


300.513 


2434.3 






622 




371 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


315 


185.125 


505.925 






661 




29 


5 


[0,5,0,0,0,0,0,0,0,0,0,0,0,0,0] 


97 


3.899 


258.275 






106 



18 
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50 


5 


[5,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


99 


4.911 


533.775 






116 




95 


6 


[3,2,0,0,0,1,0,0,0,0,0,0,0,0,0] 


115 


9.017 


73.725 






154 




96 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


120 


7.652 


24.55 






171 




99 


9 


[9,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


127 


9.08 


25.7 






170 




109 


9 


[9,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


138 


11.984 


30.225 






187 




110 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


124 


7.249 


31.025 






192 


S15(5) 


119 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


127 


9.411 


36.375 


685.921 


45.72807 


211 




131 


10 


[9,1,0,0,0,0,0,0,0,0,0,0,0,0,0] 


121 


10.939 


45.3 






252 




134 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


174 


17.188 


48.075 






258 




135 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


173 


20.678 


225.525 






288 




152 


10 


[8,2,0,0,0,0,0,0,0,0,0,0,0,0,0] 


180 


28.469 


144.375 






338 




162 


9 


[8,0,0,0,1,0,0,0,0,0,0,0,0,0,0] 


205 


32.88 


83.55 






357 




180 


10 


[8,2,0,0,0,0,0,0,0,0,0,0,0,0,0] 


269 


71.89 


269.85 






471 




201 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


306 


440.674 


1883.525 






584 




131 


8 


[0,8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


188 


35.321 


321.325 






264 




145 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


191 


34.252 


332.3 






265 




249 


9 


[6,2,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


233 


54.57 


550.725 






340 




257 


9 


[6,2,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


233 


53.913 


569.15 






352 




260 


8 


[8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


197 


47.428 


573.775 






355 




319 


9 


[1,7,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


244 


90.809 


785.925 






451 




354 


9 


[4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


256 


97.12 


938.925 






500 




459 


10 


[0,10,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


398 


182.085 


1700.525 






787 




465 


9 


[9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


356 


363.725 


1747.575 






752 


S20(l) 


469 


11 


[3,8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


317 


239.548 


1802.75 


19603.67 


980.1835 


796 




487 


9 


[9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


384 


408.842 


2011.8 






865 




572 


12 


[6,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


399 


826.33 


3290.4 






1160 




575 


10 


[0,10,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


542 


3123.4 


3414.425 






1235 




587 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


487 


1356.94 


3657.525 






1273 




606 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


507 


2100.84 


4119.5 






1389 




607 


10 


[5,5,1,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0] 


787 


2894.21 


4181.7 






1387 




652 


11 


[9,0,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


683 


2013.25 


5538.2 






1673 




674 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


821 


1999.747 


6288.875 






1732 




694 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


726 


2991.08 


7185.075 






1851 




762 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


118 


690.26 


11142.2 






2375 




57 


4 


[0,2,0,0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0] 


146 


15.003 


108.725 






186 




105 


7 


[7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


169 


20.95 


158.525 






211 




182 


9 


[6,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


195 


31.558 


311.125 






304 




186 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


272 


61.55 


328.45 






364 




201 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


319 


300.1 


374.2 






409 




204 


9 


[9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


284 


63.721 


394.55 






427 




254 


9 


[8,0,0,0,1,0,0,0,0,0,0,0,0,0,0] 


378 


153.266 


635.725 






599 




259 


10 


[7,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


295 


72.789 


653.9 






530 




263 


10 


[6,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


294 


69.81 


695.8 






612 



THE OMEGA INVARIANT BY OPTIMIZATION OVER AN EFF. SET 



19 



S20(2) 


274 


9 


[4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


336 


172.839 


751.05 


11703.9 


585.1952 


587 




275 


11 


[8,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


414 


217.261 


798.075 






702 




294 


8 


[5,1,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


359 


396.194 


915.575 






612 




295 


11 


[8,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


488 


2342.7 


975.65 






802 




298 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


463 


174.03 


990.675 






773 




307 


10 


[6,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


391 


187.431 


1084.4 






808 




338 


13 


[13,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


484 


432.837 


1546.65 






1001 




367 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


608 


2000.37 


2113.475 






1248 




393 


12 


[3,9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


658 


2531.045 


2889.325 






1502 




417 


11 


[5,0,2,0,0,0,0,0,0,2,3,0,0,0,0,0,0,0,0,0] 


563 


1093.34 


3737.325 






1686 




431 


14 


[6,8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


781 


1367.11 


* 






+ 




85 


4 


[0,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


147 


25.646 


341.175 






230 




298 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


316 


99.513 


864.5 






455 




333 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


322 


91.273 


1007.75 






465 




342 


16 


[16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


326 


99.809 


1026.075 






466 




349 


16 


[16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


307 


76.393 


1075.375 






512 




358 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


324 


86.003 


1092.975 






480 




401 


12 


[10,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


394 


154.683 


1372.975 






631 




415 


16 


[16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


448 


293.327 


1474 






687 




462 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


361 


135.922 


1827.75 






691 


S20(3) 


480 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


527 


261.581 


1982.85 


8912.075 


445.6038 


786 




556 


16 


[16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


610 


592.666 


2975.05 






1028 




569 


18 


[18,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


695 


2453.98 


3284.75 






1158 




583 


19 


[19,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


711 


1496.19 


3440.55 






1164 




609 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


57 


11.671 


4037.25 






1290 




619 


13 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


912 


990.061 


4518.725 






1386 




708 


18 


[18,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


582 


569.397 


+ 






H= 




710 


15 


[11,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


733 


673.064 


+ 






* 




752 


18 


[18,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


777 


722.174 


+ 






* 




821 


21 


[18,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0] 


108 


35.333 


+ 






* 




853 


21 


[19,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


108 


43.389 


+ 






* 




81 


5 


[0,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


140 


17.706 


219.175 






214 




107 


6 


[6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


176 


27.955 


264.375 






251 




168 


9 


[7,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


185 


26.037 


416.725 






291 




194 


9 


[6,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


239 


57.189 


527.75 






427 




230 


8 


[4,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


186 


39 


707.575 






471 




236 


9 


[7,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


198 


49.312 


735.875 






474 




274 


9 


[8,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


284 


524.051 


1027.225 






590 




277 


9 


[7,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


276 


113.266 


1079.7 






679 




286 


9 


[6,2,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


321 


676.315 


1143.675 






698 


S20(4) 


290 


8 


[1,7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


312 


1775.88 


1177.15 


11215.3 


560.7652 


630 




305 


10 


[7,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


256 


188.763 


1345.125 






683 




310 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


297 


85.818 


1407.6 






704 



20 
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348 


10 


[7,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


403 


392.226 


2039.675 






953 




351 


11 


[11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


432 


193.712 


2079.4 






949 




366 


10 


[9,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


383 


295.208 


2346.175 






912 




379 


10 


[3,7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


296 


1007.8 


2735.2 






1116 




396 


11 


[10,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


560 


3095.66 


3283.85 






1222 




416 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


541 


693.549 


+ 






* 




521 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


611 


955.771 


+ 






* 




583 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


982 


1000.085 


+ 






* 




101 


8 


[0,8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


57 


8.393 


488.75 






246 




141 


7 


[7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


146 


33.759 


567.975 






260 




279 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


99 


23.478 


1170.65 






502 




314 


10 


[10,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


199 


68.421 


1328.55 






457 




329 


11 


[7,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


85 


17.706 


1428.15 






461 




369 


11 


[5,5,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


106 


25.178 


1747.875 






493 




399 


11 


[7,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


156 


90.957 


2166.55 






711 




425 


11 


[5,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


65 


26.208 


2477.4 






718 




438 


13 


[13,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


115 


37.799 


2648.8 






732 


S20(5) 


447 


15 


[15,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


60 


11.342 


2771.675 


10007.05 


500.3524 


808 




477 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


138 


51.901 


3357.7 






929 




501 


16 


[16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


40 


7.425 


3884.325 






1026 




534 


12 


[12,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


180 


145.075 


4557.05 






983 




536 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


83 


23.15 


4752.25 






1090 




555 


13 


[9,3,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


166 


63.804 


5404.225 






1190 




574 


13 


[13,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


345 


VYY.094 


* 






* 




620 


14 


[14,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


721 


654.063 


+ 






* 




727 


18 


[18,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


882 


2140.435 


* 






+ 




786 


17 


[17,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


734 


3000.11 


+ 






* 




871 


17 


[17,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 


813 


2800.75 


* 






+ 



Table 3: Results of computational experiments for p = 15, 20. 



From the computational tests we can assert that our method improves the current algorithm imple- 
mented in GAP used to compute the w invariant of a numerical semigroup because several reasons: 

• The CPU times of our algorithm are much smaller than the brute force implementation in GAP. 
Observe that the average quotient between timej and GAPtime is 0.23. 

• Our algorithm avoids the complete enumeration of the set of non-dominated solutions of (SMIPj). 
Note that the average percentage of explored solutions (itj) with respect to the number of non- 
dominated solutions (#min) is 58.9% while in GAP the 100% of the solutions are explored. The 
algorithm proposed in [5] also needs to compute a large number of solutions to obtain the uj 
invariant, so although we have not explicitly compared the running times with this method, our 
algorithm only searches a part of the efficient solutions and then, it must take much smaller CPU 
times. 
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• Algorithm [1] is able to compute the uj invariant of the whole battery of experiments while the GAP 
method was not able to solve 14 problems (*). 

5. Conclusions 

In this paper we propose a new methodology, based on discrete optimization tools, to compute a 
specific arithmetic invariant for numerical semigroups. We use a result in |S] f Theorem 12. 2|) to translate 
the algebraic problem to the problem of optimizing a linear function over the non-dominated solutions of 
a multiobjective integer linear program that is basically a knapsack-type problem. Then, an adaptation of 
the methodology in |20| is applied to compute more efficiently than the current methods, the w invariant 
of a numerical semigroup. The key here is to avoid the complete enumeration of the non-dominated 
solutions of the multiobjective problem to solve the problem. 

The algorithm may be improved by using better bounds or by adding more information to the problem. 
It is here where bidirectional improvements in both fields, optimization and commutative algebra, may 
help to have better computational results. 

Note that the proposed methodology can be adapted conveniently to compute the w invariant of some 
other commutative semigroups. In particular, in the case when the semigroup is defined by a system of 
linear Diophantine equations, like in full and affine semigroups (also called saturated semigroups), from 
Theorem 12.21 a similar result than the one in Theorem 12.31 can be derived but where the feasible region 
consists of several equations instead of the knapsack type equation. For those semigroups, can also be 
adapted Algorithm [1] to compute the uj invariant. 

Furthermore, some other arithmetic invariant, as the tame degree, should also be modeled as a discrete 
optimization problem, and optimization techniques can help to solve more efficiently these problems. The 
analysis of different arithmetic invariant by discrete optimization is left for further research. Different 
approaches but also using tools from operational research may be used to solve problems for numerical 
semigroups such as the computation of minimal decomposition of a numerical semigroup into irreducible 
numerical semigroups or to define new indices that can be defined as optimal values of a linear objective 
function over the Kunz's polytopc. 

We hope that this paper is the first of many bidirectional collaborations between the fields of algebra 
and operation research. The two fields seem to be very different from each other, but that may help to 
get new results in both areas and to understand with different viewpoints notions that usually live in one 
of them. 
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